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Single and two-electron systems confined in nanowire quantum dots are studied in the context 
of spin-orbit coupling effects. Anisotropy of the spin-orbit interaction is discussed in terms of the 
system geometry and orientation of the external magnetic field. We find that in the limit of strong 
lateral confinement the electron spin becomes well defined in spite of the presence of spin-orbit 
coupling. We present an analytical solution for the one-dimensional limit and study its applicability 
for finite thickness nanowires by comparing the results with a full three-dimensional calculation. The 
present results are confronted with the recent measurements of the effective Lande factor anisotropy 
in InSb nanowire quantum dots [S. Nadj-Perge et al, Phys. Rev. Lett. 108, 166801 (2012)]. 
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I. INTRODUCTION 

There is a growing interest in gated semiconductor 
nanowires in the context of possible applications for spin- 
operating devices^ - — These structures provide a good ba- 
sis for creation of small quantum dots induced by external 
voltages. Energy spectra of such dots as determined 5 by 
transport spectroscopy bear distinct signatures of strong 
spin-orbit (SO) interaction which results from the struc- 
ture inversion asymmetry (Rashba SO coupling^) or the 
bulk inversion asymmetry (Dresselhaus SO interaction 7 -). 
SO coupling mixes spin and orbital degrees of freedom 
and opens the possibility of fully electrical control of the 
electron spin^i — allow for electron spin relaxation me- 
diated by phononS ) 10 ' 11 and introduce anisotropic correc- 
tions to spin SWAP process^ 

The SO coupling opens avoided crossings 5 - in the quan- 
tum dot energy spectra as function of the external mag- 
netic field (B). The width of the avoided crossings 
between energy levels of different spin states depends 
on the orientation of B vector. The dependence re- 
veals the anisotropy of the SO interaction4i^r— More- 
over, the mixing of the spin states by SO coupling de- 
termines an effective Lande factor (g-factor) and its 
anisotropy^ in function of the magnetic field orienta- 
tion. In a nanowire quantum dots the effective g-factor 
was recently measured^ in electric dipole spin resonance 
(EDSR) experiments for a two-electron spin-blocked con- 
figuration. The anisotropy of SO interaction is a relevant 
issue for spin qubit manipulation^ as well as for helical 
spin liquids^ 7 - which in the proximity of a superconductor 
can be used for observation of Majorana fcrmions 

It is well known, that in the presence of SO coupling 
the electron spin can be well defined in the stationary 
eigenstates only for equal Rashba and Dresselhaus SO 
coupling constants^ This fact was exploited in a pro- 
posal of nonballistic spin field effect transistor— and for 
prediction^ of persistent spin helixj^i In the present 
work we demonstrate that in the limit of strong lateral 



confinement the electron spins become well defined in 
the direction perpendicular to the wire axis and the ex- 
ternal electric field vector in spite of the presence of the 
pure Rashba coupling. In a general case, we show that 
depending of the external magnetic field orientation the 
electron spin states can be strongly mixed or an almost 
complete spin polarization is obtained. 

For a description of narrow nanowires a one- 
dimensional model is commonly used^ In the present 
work we present an analytical form of eigenstates for this 
approximation. The analytical form of the SO-couplcd 
wave functions explains a different strengths of the spin- 
splittings for varied orientation of the magnetic field. 
We study applicability of the one-dimensional model for 
a finite radius nanowire by comparing its results with 
three-dimensional calculation for various geometries of 
nanowire quantum dot. We relate our results in the two- 
electron regime to the measurements of effective g-factor 
on an InSb nanowire quantum dot in the EDSR experi- 
ment of Ref. H. 



II. THEORY 

We consider single-electron nanowire quantum dot de- 
scribed by the three-dimensional Hamiltonian, 



H 



2m* 



V(r) + H SO + ^9^bB ■ a 



lelF ■ 



(1) 



where k = — iV + eA/h with the gauge A = 
B(z sm<p, 0, ycos </>). The magnetic field is aligned in the 
xy-plane with an angle </> between the B and a;-axis - 
in such a case the Zeeman term stands i<?/i£B • a = 
\[iBgB(a x coscf) + (j y sin0). We account for Rashba SO 
coupling Hso = a o ■ (c X k) as the main SO in- 
teraction type in the [111] grown InSb nanowires^ V(r) 
stands for the confinement potential which we take in a 
separable form V(r) = Vi(y, z)V(x) where Vi(y, z) is a 
400 meV deep two-dimensional circular quantum well of 
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radius R, and V(x) is a infinite quantum well with width 
L. We assume electric field F = (0,0,^) with non-zero 
component in the z-direction (perpendicular to the axis 
of the wire) due to the gating of the nanowire>i~— For this 
form of the electric field if so = ct(a x k y — cr y k x ) where 
a = a F z . 

We solve the Schrodinger equation using finite dif- 
ference method on a grid employing method previously 
used for self-organized quantum dots^ where the single- 
electron solutions are found in basis constructed from 
one-dimensional spin-orbitals. 

The two-electron system is described by the Hamilto- 
nian, 



H = hi 



4irs e\ri - r 2 | 



(2) 



and calculation is performed within configuration inter- 
action scheme which treats the Coulomb interaction in 
exact manner including electron-electron exchange and 
correlation. 

We adopt material parameters 2 ^ for InSb, namely : 
to* = 0.0147«o, g = —51, e = 16.5, and ao = 5 nm 2 . In 
the bulk of the paper we choose F z = 50 kV/cm which 
results in SO interaction strength a = 25 meVnm. Unless 
stated otherwise we take L = 300 nm. 



III. RESULTS 

A. Single electron in a finite thickness nanowire 
quantum dot 

Lowest part of the energy spectrum of the single- 
electron quantum dot is presented in Fig. [TJ In the ab- 
sence of the magnetic field all the levels are Kramer's dou- 
blets. We include residual magnetic field with B = 5 mT 
and check the spin polarization along the magnetic field 
direction (calculated as (sb) = (s x ) cos(0) -I- (s y ) sin(0)). 
In Fig. [2ja) we present that the spin polarization un- 
dergoes oscillatory changes as function of B orientation. 
This reflects the presence of the easy and hard spin po- 
larization axes in the system. For the magnetic field ori- 
ented perpendicular to the nanowire axis the spin is easily 
polarized - taking values close to 1 [fi/2]. On the other 
hand for B oriented along the wire the (sb) is around 
0.885 h/2. The amplitude of the oscillations depends on 
the nanowire radius (compare the curves in Fig. [5{a) for 
three values of R) and the oscillations are the strongest 
for narrow nanowire with R = 10 nm. Let us set <j) = 90° 
and further inspect the spin polarization. In Fig. [2jb) we 
plot mean value of the spin-y component in function of 
the wire radius R. We observe that as the wire becomes 
narrower the spin polarization becomes almost complete 
(i.e., 1 — {s y )2/% < 10~ 4 for R = 1 nm) despite the pres- 
ence of the SO coupling. On the other hand as the wire 
becomes wider the spin-polarization drops with the slope 
of the curves in Fig. [5] depending on the SO strength a. 
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FIG. 1. (color online) Single-electron energy spectrum for the 
SO coupled nanowire quantum dot with radius R = 50 nm 
and SO interaction strength a = 25 meVnm plotted with lines 
for two orientations of the magnetic field. The crosses are the 
results obtained from asymptotic one-dimensional solution - 
see text. With | t) an d | 4-) we mark the spin-polarization 
of the states parallel and antiparallel to the magnetic field 
respectively as found without SO coupling. 



When the magnetic field is increased it splits the dou- 
blets - see the energy levels in Fig. [1] The energy split- 
tings obtained for magnetic field perpendicular to the 
nanowire axis (red curves in Fig. [T]) are stronger than the 
ones obtained for magnetic field parallel to the nanowire 
axis (black curves in Fig. [IJ. In the following we explain 
this observation. 



B. Asymptotic solution (ID limit) 

When the wire becomes narrow the energy of the states 
excited in the radial direction (in our case the y and z 
directions) rises. It is reasonable then to inspect the case 
where the radial degrees of freedom are decoupled from 
the longitudinal one (the x-direction) . Such a system is 
described by the one-dimensional Hamiltonian ) 22 ! 23 
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hiD = 



2m* 



V{x)-aa y k x + \-n B gB{a x cos (j>+a y sine 



( 3 ) 

where k x = —i-j^.- 

Generally, the analytical solution for a SO coupled con- 
fined systems are not known with the exception of a spe- 
cial case of equal strength of Rashba and Dressclhaus 
coupling described in Rcf. Il9l . Here we note however 
that in the absence of the magnetic field {B = 0) the 
Hamiltonian ([3]) commutes with spin-y Pauli matrix and 
its eigenstates are the states with well defined spin in 
the y-direction. We find that for a quasi one-dimensional 
nanowire the spin-orbitals (where N stands for the or- 
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FIG. 2. (color online) (a) Mean value of the spin along mag- 
netic field direction obtained for the nanowire with radius 
R = 10 nm (blue curve), R — 50 nm (green curve) and 
R = 100 nm (black curve), (a) mean value of the spin-i/ com- 
ponent for magnetic field aligned along y-direction in function 
of the nanowire radius R. (a) and (b) are obtained for B — 5 
mT. Results for a = 50 meVnm are obtained with F- = 50 
kV/cm. 



bital quantum number and ± for the spin polarization of 
the state) have the form, 



1 



N± 



V2\± 
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<^Ar(a;)exp 



lam 
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(4) 



where ip^{x) are spin-independent cigenstates of Hamil- 
tonian ([3]) for a = and B — 0. The eigenenergies of 
the Hamiltonian ([3]) are E\d = -Ea=o.N + Eso where 
Eso = —a 2 m*/(2h 2 ) is the energy shift to the whole 
energy spectrum introduced by the SO interaction^ and 
E a= o t N is an energy level of the N'th eigenstatc obtained 
without SO coupling. 

The magnetic field affects the energy levels of a 
strongly confined electron mainly through the Zeeman 
spin-splitting. To investigate its influence on the SO 
cigenstates with orbital excitation TV let us diagonalize 
hio for B > in a basis consisting of a degenerate pair 
v E , tv+ and Vl/jv-- The Hamiltonian matrix is 



(y N+ \h 1D \^ N+ ) (y N -\h 1D \^ N+ ) 
(t? N+ \h 1D \y N -) (y N J\h 1D \V N -) 



where the diagonal elements are defined as follows 



(^ N± \hi D \^ N± ) = E XD ± -g(i B Bsm<f>, 



(5) 
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dx cos < 
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Let us denote = 
iJ\(p N \ 2 sm(^^x)dx 



- J \Vn\ 2 cos( 2 'p l x)dx and n N 



The eigenstates of the matrix (J5j arc 



E 



N± 



Eid ± 



\ 2 



N 



)cos 2 0- (8) 



The energy difference between the states depends on 
the orientation of the magnetic field (angle <fi) as well 
as the parameters Xn and kn that control strength of 
the anisotropy of the spin splittings for rotated magnetic 
field. For the symmetric infinite quantum well confine- 
ment along the wire (cc-direction) we obtain^ 



Ai 



h b Tr 2 sa\{Lam*/h z ) 
am* L(tt 2 Ti a — a 2 m* 2 L 2 ) ' 



and 



4h 6 n 2 sin{Lam*/Ti 2 ) 



a 2 m* 2 



L 2 ) 



(9) 



(10) 



and Ki = K2 = for the two lowest orbital excited states. 
The Ajv depends on the nanowire length and the SO 
strength. In Fig. [3] we present the Ai parameter in func- 
tion of L and a. With the light-green dashed curve we 
depict the SO length l so = h/ (m*a). We observe that Ai 
drops quickly when the length of the dot becomes greater 
than SO length. The shape of the Ai dependence on the 
SO strength for different quantum dot lengths is pre- 
sented in Fig. [3jb) showing that the SO effects strongly 
depends on the quantum dot geometry and that Ai goes 
to 1 for vanishing SO coupling 

The smaller Atv is the stronger effect of SO coupling 
are. In particular for the magnetic field parallel to 
the nanowire axis the energy of the spin splitting is 
Es = gfisBX^- Consequently the splitting can even 
go to due to strong mixing of the spin states by the SO 
interaction (the light blue region in Fig. EJa)). 

When the magnetic field is aligned in the direction per- 
pendicular to the nanowire axis, i.e. = 90° or = 270° 
the off-diagonal elements of matrix ([S]) vanish and the 
energy levels are split by Zeeman energy with the bulk 
value of the g-factor. This is the reason for stronger spin 
splittings of the red curves in Fig. [TJ For this configura- 
tion the spin-orbitals are separable into spin and orbital 
parts despite the presence of SO interaction and they 
have exact form of Eq. (£[]). For any other orientation 
of the magnetic field the off-diagonal elements mix the 
cigenstates This results in decreasing the spin split- 
tings by the SO interaction by the amount that depends 
on the Ajv and kn parameters - the spatial extent of the 
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FIG. 3. (color online) Parameter Ai in function of dot length 
L and SO coupling strength a. (b) cross-section of (a) for 
three different dot lengths L. 



wave function along the nanowire and the strength of 
the SO coupling. Moreover the electron spin is no longer 
well defined as the wave function is not separable into 
spin and orbital degree of freedom. 

We plot the energy spectrum obtained from Eq. © 
(shifted to match the energies obtained in the three- 
dimensional calculation at B — 0) with the crosses on 
Fig. [TJ The spin splitting obtained from the one- 
dimensional model well describes the results of the three- 
dimensional calculation. The only discrepancy is visible 
for the energy levels of the first and the second excited 
states for B > 100 mT which is due to mixing of this two 
states by the SO interaction. 
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FIG. 4. (color online) Effective g-factor obtained for a 
nanowire with R — 10 nm (red curve), R = 50 nm (blue 
curve) and R = 100 nm (green curve) and from the one- 
dimensional model (black crosses). Results are obtained for 
B = 100 mT. 



As the magnetic field is rotated between the easy and 
hard axes the spin-polarization of the states changes 
which results in changes of the spin-splitting strength. 
The latter term in Eq. ([5]) introduces Zeeman energy 
splitting between the energy levels of the two states. We 
can see that 



g* = . 9v /1-(1-A 2 + k 2 )cos 2 



(11) 



is an effective g-factor that is dependent on the orienta- 
tion of magnetic field with the angle <j>. With crosses in 
Fig. S]we plot effective g-factor as obtained from Eq. (fTTj) 
along with the values obtained in the three-dimensional 
calculation (calculated as g* = AE/hbB, where AE is 
the energy difference between the energy of first excited 
state and the ground-state - see Fig. [T]) for different 
nanowire radii. For the nanowire radius R = 10 nm they 
match perfectly. For the other two values of the R the 
shape of the dependence holds, only the amplitude is dif- 
ferent, with the biggest discrepancy for the wide nanowire 
with R = 100 nm. 



C. Two-electron results 

The anisotropy of g-factor in the two-electron regime 
is extracted from the slopes of resonance lines in EDSR 
experiments on double quantum dotsi^ Figure [Ha) 
presents two-electron energy spectrum of weakly cou- 
pled quantum dots in a nanowire obtained in the three- 
dimensional calculation for the nanowire with radius 
R = 30 nm. Results for the magnetic field oriented 
along the nanowire axis with = (perpendicular to 
the nanowire with <f> = 90°) arc plotted with solid (dot- 
ted) curves. The confinement potential includes now a 
potential barrier of 60 nm width that separates the elec- 
trons in adjacent dots both of 120 nm width. At B = 
the ground-state is a singlet state (| tl) — I It)) split 
from the degenerate triplet states (see the inset to Fig. 
G3a)). We tune the barrier height to 5 meV to match 
the singlet-triplet separation of ~ 5 ^eV as measured 
in Rcf. y|. At B = 3 mT an avoided crossing between 
the two lowest energy levels appears for = due to 
spin mixing by the SO interaction. After the anticross- 
ing the magnetic field splits the energy levels of the two 
spin-polarized triplet states (| tt) and | 44-)) by the Zee- 
man energy. The blue and red curves in Fig. [5ja) whose 
energy does not change (after the anticrossing) with B 
are the singlet (| tl) - | It)) and triplet (| f4) + I It)) 
states with zero spin-component in the direction along 
the magnetic field. Those levels are split by exchange 
interation^ (additional splitting to those two energy lev- 
els occurs when when the g-factor along the structure is 
not constant^— ). 

The change of the magnetic field orientation (angle (j>) 
results in: i) change in strength of the spin polarization 
of the triplet states | tt) > I 11) which results in a change 
of the slope of the corresponding energy levels and ii) 
change in the exchange energy (spacing between energy 
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FIG. 5. (color online) (a) Two-electron energy spectrum of 
a coupled nanowire quantum dots with radius R — 30 nm. 
Solid curves present results for cj> = and dotted curves for 
= 90°. With | tt), 1 44), I 4-t), I U) - I 4-t) and 1 14-> + I 4-t) 
we mark the spin configuration of the states parallel or an- 
tiparallel to the magnetic field as found without SO coupling. 
The inset presents the energy levels for low values of the mag- 
netic field where the avoided crossing appears, (b) with the 
curves - effective g-factor calculated from the energy split- 
tings between the ground-state energy level and the energy 
levels depicted with blue (AEi) and red (AE2) curves in (a) 
for B — 200 mT. The circles correspond to the effective one- 
electron g-factor as obtained from Eq. [TT] (shifted down by 1) 
for a single quantum dot with length L — 120 nm. 



levels of 1 14-) - I It) an d | tl) + 1 It) states, plotted with 
blue and red curves in Fig. [5{a)). These two effects lead 
to a change of the effective g-factor which we calculate 
from the energy splittings between the ground-state and 
the first and second excited states and plot in Fig. [SJb) 
with the blue and red curves respectively. We find that 
the shape of both curves in Fig. [SJb) match the shape of 
the single-electron dependence presented in Fig. 2] only 
the amplitude of the oscillations is lower. As described 
by Eq. (TIT]) for the single-electron case the amplitude 
of g-factor oscillations depends on the dot length. In 
the present case each of the coupled quantum dots has 
a length of L = 120 nm. The effective g-factor obtained 
for a single dot of this length as calculated from Eq. (TiT]) 
is plotted with circles in Fig. [SJb). Obtained oscillations 
have similar amplitude to the ones obtained for the two- 
electron system. This suggest that the low amplitude 



in the two-electron case results from the fact that each 
electron resides in a separate dot and the shape of the os- 
cillations is controlled mainly by the single-electron spin 
polarization anisotropy. 

The almost complete spin polarization occurring for 
the magnetic field oriented perpendicular to the nanowire 
results in closing of the the avoided crossing between the 
energy levels in the energy spectrum (see dotted curves in 
Fig. HJa)). This finding is in agreement with the experi- 
mental results [||. Also the shape of the g-factor depen- 
dence is similar to the one obtained in the experiment. 
In particular an agreement is obtained in the context of 
the slight change of the oscillation amplitude of the red 
and blue curves in Fig. [5jb). This difference in ampli- 
tudes is due to modification of the exchange energy that 
separates the energy levels of the singlet (| fj,) — I tt)) 
and triplet (| f4-) + 1 tt)) states by rotated magnetic field. 
However, the experimental dependence of the effective g- 
factor is shifted (with minima at <fi — 124° and <j> = 304°) 
with respect to the present result. We performed calcu- 
lations for quantum dots in a nanowire of larger radius 
(R = 100 nm) ruling out the possible orbital effects of 
the magnetic field as a reason for the shift. Our previ- 
ous studyi^ indicates that any additional electric fields 
in the device activating additional Rashba terms, or the 
presence of Dresselhaus coupling would change the mag- 
netic field angle for which the avoided crossing vanishes^ 
along with the shift of the effective g-factor angular de- 
pendence. 



IV. SUMMARY AND CONCLUSIONS 

In the present work we studied the anisotropy of the 
spin polarization in a narrow nanowire quantum dot in 
the presence of SO coupling. By the solution of the 
three-dimensional Schrodingcr equation we showed that 
strength of the spin polarization under the presence of 
Rashba SO interaction depends on the orientation of the 
magnetic field and that there are hard and easy spin po- 
larization axes as SO interaction tends to polarize spins 
in the direction perpendicular to the nanowire. For the 
magnetic field aligned in this direction the electron spin 
polarization can be nearly complete depending on the 
nanowire radius. We presented analytical solution for 
the one-dimensional limit and compared its results with 
the calculation for a finite thickness nanowire. The spin 
polarization anisotropy results in the effective g-factor 
dependence on the magnetic field orientation. The form 
of this dependence matches the one obtained in the ex- 
periment. 
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